Gate-dependent spin-orbit coupling in multi-electron carbon nanotubes 
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Understanding how the orbital motion of electrons is coupled to the spin degree of freedom in 
nanoscale systems is central for applications in spin-based electronics and quantum computation. 
We demonstrate this coupling of spin and orbit in a carbon nanotube quantum dot in the general 
multi-electron regime in presence of finite disorder. Further, we find a strong systematic depen- 
dence of the spin-orbit coupling on the electron occupation of the quantum dot. This dependence, 
which even includes a sign change is not demonstrated in any other system and follows from the 
curvature-induced spin-orbit split Dirac-spectrum of the underlying graphene lattice. Our findings 
unambiguously show that the spin-orbit coupling is a general property of nanotube quantum dots 
which provide a unique platform for the study of spin-orbit effects and their applications. 
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The interaction of the spin of electrons with their orbital 
motion has become a focus of attention in quantum dot 
research. On the one hand, this spin-orbit interaction 
(SOI) provides a route for spin decoherence, which is 
unwanted for purposes of quantum computation [H-Q- 
On the other hand, if properly controlled, the SOI can 
be utilized as a means of electrically manipulating the 
spin degree of freedom 0-0|- 

In this context, carbon nanotubes provide a number 
of attractive features, including large confinement en- 
ergies, nearly nuclear-spin-free environment, and, most 
importantly, the details of the energy level structure is 
theoretically well understood and modeled, as well as 
experimentally highly reproducible. Remarkably, the 
SOI in nanotubes was largely overlooked in the first 
two decades of nanotube research and was only recently 
demonstrated by Kuemmeth et al. for the special 
case of a single carrier in ultra-clean CNT quantum 
dots@, H, 0- Except for these reports, the SOI in 
nanotubes is experimentally unexplored. Theoretically, 
the focus has exclusively been on the SOI-modified band 
structure of disorder-free nanotubes [l0rll3 | . Therefore, 
two important questions remain: how the effective SOI 
depends on electron filling and how it appears in the 
general case of quantum dots subject to disorder. Here 
we answer these two questions. 

Firstly, by low-temperature electron transport we 
demonstrate the presence of a significant SOI in a 
disordered CNT quantum dot holding hundreds of 
electrons. We identify and analyze the role of SOI in 
the energy spectrum for one, two, and three electrons 
in the four-fold degenerate CNT electronic shells, thus 
describing shells at any electron filling. By rotating the 
sample, we present for the first time spectroscopy of the 
same charge-states for magnetic fields both parallel and 
perpendicular to the nanotube axis, thus controlling the 
coupling to the orbital magnetic moment. Remarkably, 



a single-electron model taking into account both SOI 
and disorder quantitatively describes all essential details 
of the multi-electron quantum dot spectra. Secondly, 
by changing the electron occupancy we are able to tune 
the effective SOI and even reverse its sign in accordance 
with the expected curvature-induced spin-orbit splitting 
of the underlying graphene Dirac spectrum [l], IllHlij . 
Such systematic dependence has not been demon- 
strated in any other material system and may enable 
a new range of spin-orbit related applications. This 
microscopic understanding and detailed modeling is in 
stark contrasts to situations encountered in alternative 
strong-SOI quantum-dot materials, such as InAs or 
InSb nanowires, where the effective SOI arises from bulk 
crystal effects combined with unknown contributions 
from surface effects, strain and crystal defects (l5j. 
These systems often exhibit semi-random fluctuations 
of, e.g., the (/-factor as single electrons are added(l6lll7j. 
Thus, beyond fundamental interest and the prospect 
of realizing recent proposals of SOI-induced spin con- 
trol in CNTsp], [l8j], our findings pave the way for new 
designs of experiments utilizing the SOI in quantum dots. 



Zero-field splitting of the four-fold degeneracy 

Our experimental setup is presented in Fig. la. We 
fabricate devices of single-wall CNTs on highly doped 
Si substrates capped with an insulating layer of Si02 
(see Methods section). The size of the quantum dots 
is defined by the contact separation (400 nm) and the 
electrical properties are investigated in a voltage biased 
two-terminal configuration applying a voltage V s d be- 
tween source-drain contacts and measuring the resulting 
current /. The differential conductance dI/dV s d is mea- 
sured by standard lock-in techniques. When biased with 
a voltage V g , the Si substrate acts as an electrostatic 
gate controlling the electron occupancy of the dot. The 
devices are measured at T = 100 mK in a 3 He/ 4 He 



2 



dilution refrigerator, fitted with a piezo-rotator allowing 
in-plane rotations of the device in magnetic fields up to 
9T. 

Figure lb shows a typical measurement of dI/dV sc i 
vs. V sc i and V g in the multi-electron regime of a small- 
band-gap semiconducting nanotube. The pattern of 
diamond-shaped regions of low conductance is expected 
for a quantum dot in the Coulomb blockade regime 
and within each diamond the quantum dot hosts a 
fixed number of electrons N, increasing one-by-one with 
increasing V g . The energy -E a dd, required for adding a 
single electron appears as the diamond heights and has 
been extracted in Fig. lc. The four-electron periodicity 
clearly observed in Fig. lb and c reflects the near four- 
fold degeneracy in the nanotube energy spectrum fl9l I20I] : 
one factor of 2 from the intrinsic spin (t,i) and one 
factor of 2 from the so-called isospin (K, K') that stems 
from the rotational symmetry of the nanotube - electrons 
orbit the CNT in a clockwise or anticlockwise direction. 
As is generally observed [l9l - l23| . the addition energy for 
the second electron in each quartet (yellow in Fig. lc) 
exceeds those for one and three. This was previously 
interpreted as a result of disorder-induced coupling 
Akk 1 of the clockwise and anticlockwise states [2l], [24[ 
that splits the spectrum into two spin-degenerate pairs 
of bonding/antibonding states separated by Akk 1 - As 
mentioned, Kuemmeth et al. showed recently that for 
the first electron in an ultra-clean suspended nanotube 
quantum dot, the splitting was instead dominated by the 
spin-orbit coupling. The first question we address here 
is whether SOI also appears in the many-electron regime 
and how it may be modified or masked by disorder. 



Modeling spin-orbit coupling and disorder 



zero-field spectrum is also split into two doublets, but 
the field dependence is markedly different and no avoided 
crossing appears. In the simplest picture, this behavior 
originates from coupling of the electron spin to an 
effective magnetic field B50 = — (v x E)/c 2 experienced 
by the electron as it moves with velocity v in an electric 
field E. Here the speed of light, c, reflects the relativistic 
origin of the effect. In nanotubes, the curvature of the 
graphene lattice generates an effective radial electric 
field, and since the velocity is mainly circumferential 
(and opposite for K and K' ), B50 polarizes the spins 
along the nanotube axis and favors parallel alignment 
of the spin and orbital angular momentum. Thus, even 
in the absence of disorder, the spectrum splits into two 
Kramers doublets (K |, K' t) and (K t, K' |) separated 
by Aso- Interestingly, since a perpendicular field does 
not couple K and K' the doublets do not split along 
B± on the figure. As a consequence, the g-factor, when 
measured in a perpendicular field, will vary from zero 
when Aso ^ Akk' (Fig. 2c) to 2 in the opposite limit 
(Fig. 2b). 

The final case, including both disorder and SOI, is 
of particular importance for the present study, and 
the calculated spectrum is displayed in Fig. 2d for 
Akk' > Aso- Importantly, the effects of SOI are not 
masked despite the dominating disorder: For parallel 
field, SOI remains responsible for an asymmetric split- 
ting of the Kramers doublets (a, /3) vs. (5, 7), and the 
appearance of an additional degeneracy in the spectrum 
at finite field (S and 7 states). In a perpendicular field, 
the effect of SOI is to suppress the Zeeman splitting of 
the two doublets and since the eigenstates of the SOI 
have spins along the nanotube axis it couples the states 
with spins polarized along B± resulting in the avoided 
crossing indicated on the figure. 



Performing level spectroscopy with a magnetic field B 
applied either parallel (Bn) or perpendicular (B±) to the 
nanotube axis proves to be a powerful tool to analyze 
the separate contributions from disorder and spin-orbit 
coupling. This is illustrated in Figs. 2a-d, which show 
calculated single-particle energy level spectra for four 
limiting combinations of Akk' and the effective spin- 
orbit coupling Aso (all limits are relevant for nanotube 
devices depending on the degree of disorder and CNT 
structure [lj, Hli details of the model are provided in the 
Supplementary Information). In all cases a parallel field 
separates the four states into pairs of increasing (X-like 
states) and decreasing (if'-like states) energies. The 
magnitude of the shift is given by the orbital g-factor 
5orb reflecting the coupling of B\\ to the orbital magnetic 
moment caused by motion around the CNT (25[. Further, 
each pair exhibits a smaller internal splitting due to the 
Zeeman effect. Figure 2b shows the disorder induced 
coupling of K and K' states resulting in an avoided 
crossing at Bh = and the zero-field splitting discussed 
above. In the opposite limit with SOI only (Fig. 2c), the 



Spin-orbit interaction revealed by spectroscopy 

With this in mind we now focus on the quartet with 
4iVo « 180 electrons highlighted in Fig. lb and expanded 
in Fig. 3a. In order to investigate the level structure we 
perform cotunneling spectroscopy, as illustrated in the 
schematic Fig. 3bj26|: In Coulomb blockade, whenever 
e I V^d I matches the energy of a transition from the ground 
state a to an excited state (/3,7,<5), inelastic cotunnel 
processes, that leave the quantum dot in the excited 
state, become available for transport. This significantly 
increases the current and gives rise to steps in the 
conductance. These appear as gate-independent features 
in Fig. 3a (arrows) and are clearly seen in the inset 
showing a trace through the center of the one-electron 
(4A^o + 1) diamond along the dashed line. Thus following 
the magnetic field dependence of this trace, as shown 
in Fig. 3c, maps out the level structure. The energies 
of the excitations are given by the inflection points of 
the curve (i.e. peaks/dips of d 2 1 / dV 2 d )]21§ and the level 
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evolution is therefore directly evident in Fig. 3d-i, which 
show color maps of the second derivative vs. V s d and 
B±,B\\ for V g positioned in the center of the one, two, 
and three electron charge states. As explained below, 
the SOI is clearly expressed in all three spectra. 

Consider first the one-electron case: In a parallel field 
(Fig. 3d), the asymmetric splitting of the two doublets 
is evident (black vs. green arrows) and applying the 
field perpendicularly (Panel g), the SOI is directly 
expressed as the avoided crossing indicated on the 
figure. The measurement is in near-perfect agreement 
with the single-particle excitation spectrum calculated 
by subtracting the energies of Fig. 3b and shown by 
the solid lines. The calculation depends on only three 
parameters: Aso = 0.15 meV set directly by the 
avoided crossing, Akk' = 0.45 meV determined from 
the zero-field splitting of the doublets (see Fig. 2d), and 
<?orb = 5.7 set by the slopes of the excitation lines from 
a to 7, S in Panel d. 

Consider now the role of SOI for the doubly occupied 
CNT quartet. This situation is of particular importance 
for quantum computation as a paradigm for preparation 
of entangled states and a fundamental part of Pauli 
blockade in double quantum dots. 2$]. Figures 3e,h show 
the measured spectra in parallel and perpendicular 
fields. The model perfectly describes the measurement 
and now contains no free parameters since these are 
fixed by one-electron measurement. Six states are 
expected: the ground state singlet-like state So formed 
by the two electrons occupying the low-energy Kramers 
doublet [29|. three triplet-like f-,T ,T+ and a singlet- 
like state Si, which all use one state from each doublet 
and the singlet-like S2 with both electrons occupying 
the high-energy doublet. The ground state §0 does 
not appear directly in the measurement, but sets the 
origin for the cotunneling excitations. The excitation 
to the high-energy §2 (dashed line) is absent in the 
experiment, since it cannot be reached by promoting 
only a single electron from So. In Fig. 3h excitations to 
T_ and T + are clearly observed, while S\ and To merge 
into a single high-intensity peak showing that any ex- 
change splitting J is below the spectroscopic line width 
~ 100^eV[23j. In other quartets an exchange splitting 
is indeed observed (see Supplementary Information). In 
the two-electron spectra the SOI is directly expressed as 
the avoided crossing at B± m 4.5 T accompanying the 
So <H> T_ ground state transitionpjl Ho[. In quartets 
of yet stronger tunnel-coupling it is replaced by a 
singlet-triplet Kondo resonance [3 1| (see Supplementary 
Information) . 

Finally, the spectrum of three electrons in the four- 
electron shell is equivalent to that of a single hole 
in a full shell; at low fields the 5-state becomes the 
ground-state and 7 the first excited state, while a and 
/3 then constitute the excited doublet. As seen by 
comparing Fig. 2b and 2d SOI breaks the intra-shell 
electron-hole symmetry of the nanotube spectrum. This 
is evident in the experiment when comparing Fig. 3d 



and 3f: In 3f, increasing B\\ , the lowest excited state 
7, barely separates from the ground state 6 and at 
Bh = 1.1 T they cross again, causing a ground state 
transition. At the crossing point, the spin-degenerate 
ground-state results in a zero-bias Kondo peak (see 
inset) (32J. Interestingly, this degeneracy also forms the 
qubit proposed in Ref. 1. For the -Bj_-dependence the 
one- and three-electron cases remain identical and Fig. 
3i exhibits again the SOI-induced avoided crossings. 



Gate-dependent spin-orbit coupling 

Having established the presence of SOI in the general 
many-electron disordered quantum dot, we now focus on 
the dependence of Aso on the quantum dot occupation. 
To this end, we have repeated the spectroscopy of Fig. 
3 for a large number of CNT quartets and in each case 
extracted Aso by fitting to the single-particle model (all 
underlying data are presented in the Supplementary In- 
formation). Figure 4a shows the result as a function of 
V g . An overall decrease of Aso is observed as electrons 
are added to the conduction band and, interestingly, a 
negative value is found in the valence band (i.e. SOI 
favoring anti-parallel rather than parallel spin and or- 
bital angular momentum, thus effectively interchanging 
the one and three-electron spectra). 

The magnitude of Aso is given by the spin-orbit split- 
ting of the underlying graphene band structure as we 
will now discuss. For flat graphene this splitting is very 
weak graphene _ x ^y) [JJ as it jg sec0 nd-order in 
the already weak atomic SOI of carbon Ag Q ~ 8meV. 
In nanotubes, however, the curvature induces a coupling 
between the 7r- and c-bands and generates a curvature- 
induced spin-orbit splitting, which is first order in the 
atomic SOI and thus greatly enhances Aso- Around a 
Dirac-point of the Brillouin zone (e.g. K) the graphene 
band structure appears as on Fig. 4b [llT - [l4l ]: The spin- 
up and spin-down Dirac cones are split by SOI both in 
energy and along k± , the momentum in the circumferen- 
tial direction of the CNT. The schematic also highlights 
the CNT band structure (Fig. 4a upper inset) obtained 
by imposing periodic boundary conditions on k±. In a 
finite-length CNT quantum dot also the wavevector along 
the nanotube axis, feii , is quantized, and letting e^v denote 
the energy of the N th longitudinal mode the effective SOI 
for a small-gap CNT becomes 

A so ,± = E« - Ef (1) 

= 2A° so ± ^(A g + A so y + e% 

T ^(A g -A so y + e%. 

Here the upper (lower) sign refers to the conduc- 
tion (valence) band, A g is the curvature induced en- 
ergy gap[33|, [34]], and the two terms A so and A so 
are the band structure spin-orbit parameters due to 
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curvature [lll - ll4l |. The separate contributions of the two 
terms are illustrated in the lower inset to Fig. 4a. Ag Q 
was found already in the work of Ando[l(| and accounts 
for the /^-separation of the Dirac-cones in Fig. 4b. For 
the CNT band structure this term acts equivalently to 
an Aharonov-Bohm flux from a parallel spin-dependent 
magnetic field, which changes the quantization conditions 
in the k± -direction. Characteristically, its contribution 
to Ago decreases with the number of electrons in the 
dot (ejv) and reverses sign for the valence band. This 
contrasts the ejv-independent contribution from the re- 
cently predicted Ag -term [l2T - [l4T ]. which acts as an effec- 
tive valley-dependent Zeeman term and accounts for the 
energy splitting of the Dirac-cones in Fig. 4b. For the 
nanotube studied in Kuemmeth et al. [9] Aso has the 
same sign for electrons and holes, i.e., \A a so \ > |A<j |. 
In our case, the negative values measured in the valence 
band demonstrate the opposite limit |Ag | < |A<j |. 
Thus, the measured gate-dependence of Aso agrees with 
the spin-orbit splitting of the graphene Dirac spectrum 
caused by the curvature of the nanotube, and fitting to 
Eq. 1 (Fig. 4a, dashed line) yields A^ c = 10±10/xeV and 
A so = 220 ± 25 ^eV[35j|. Band structure models relate 
these parameters to the structure of the nanotube (T3- 
03: A^ G = A Ag o A g D and A^ Q = \iAg /D, where 
D is the nanotube diameter and Ai ; 2 constants that de- 
pend on the CNT class (semiconducting, small-band- 
gap). Typically CVD-grown single wall CNTs have di- 
ameters in the range 1-3 nm, while obtaining D from the 
measured values of <7orb[36| gives D ss 6nm. Thus taking 
D = 1-6 nm we estimate Ai = 0.03-0.17 nm and Ao = 0.7- 
4 • 10~ 6 (nm ■ meV) -1 . While Ai is consistent with the 
value quoted in the theoretical literature (0.095 nm) [3], 
the calculated Ao value —4 • 10~ 3 (nm • meV) -1 does not 
match the experiment and similar deviations appear [l3T| 
when comparing the theory to the SOI values measured 
in Ref. 8. The origin of this discrepancy remains un- 
known, and further work on SOI in nanotubes with 
known chirality is needed to make further progress. 



Methods 

The devices are made on a highly doped Silicon wafer 
terminated by 500 nm of Si02- Alignment marks (Cr, 
70 nm) are defined by electron beam lithography prior 
to deposition of catalyst islands made of Iron nitrate 
(Fe(N03)^)., Molybdenum acetate and Alumina support 
particles [37j. The sample is then transferred to a furnace, 
where single wall carbon nanotubes are grown by chem- 
ical vapor deposition at 850-900°C in an atmosphere of 
hydrogen, argon and methane gases. Pairs of electrodes 
consisting of Au/Pd (40/10 nm) spaced by 400 nm are 
fabricated alongside the catalyst islands by standard elec- 
tron beam lithography techniques. Finally, bonding pads 
(Au/Cr 150/10nm) are made by optical lithography and 
the devices are screened by room- and low-temperature 
measurements. 

We measured the sample in an Oxford dilution refrig- 
erator fitted with an Attocube ANRv51 piezo rotator 
which allows high precision in-plane rotation of the sam- 
ple in large magnetic fields. The rotator provides resis- 
tive feedback of the actual position measured by lock-in 
techniques. For electrical filtering, room-temperature 7r- 
filters and low-temperature Thermocoax are used. The 
base temperature of the modified refrigerator is around 
100 mK. The CNT measurement setup consists of a Na- 
tional Instrument digital to analog card, custom made 
optically coupled amplifiers, a DL Instruments 1211 cur- 
rent to voltage amplifier and a Princeton Applied Re- 
search 5210 Lock-in amplifier. Standard dc and lock-in 
techniques have been used to measure current and differ- 
ential conductance dI/dV s d while d 2 I/dV^ d is obtained 
numerically 
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FIG. 1: Four-fold periodic nanotube spectrum, a, Schematic illustration of the device and setup. CNT quantum dots are 
measured at T — 100 mK in a standard two-terminal configuration in a cryostat modified to enable measurements in a high 
magnetic field at arbitrary in- plane angles 8 to the CNT axis, b, Typical measurement of the differential conductance dl /dV s d 
vs. source-drain bias V s d and gate voltage V g for a multi-electron CNT quantum dot. c, Addition energy as a function of V g . 
In b and c the characteristic filling of four-electron shells is clearly seen. 
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FIG. 2: Role of spin-orbit interaction and disorder for the nanotube energy spectrum. Calculated single-particle energy 
spectrum as a function of magnetic field applied perpendicular (B±) and parallel (-By) to the CNT axis in the limiting cases of 
neither SOI nor disorder a, disorder alone b, SOI alone c, and the two combined A KK i > Aso > d. Depending on the CNT 
type, electron filling and degree of disorder, all four situations can occur. 
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FIG. 3: Spin-orbit interaction in a disordered multi-electron nanotube quantum dot. a, Measurement of dl /dV s d vs. V s d and V g 
corresponding to the consecutive addition of four electrons to an empty shell (indicated on Fig. lb). A strong tunnel coupling 
results in significant cotunneling which is evident as horizontal lines truncating the diamonds (arrows). The black trace shows 
a cut along the dashed line, b, Schematic illustration of the relevant inelastic cotunneling processes, c, Traces along the 
dashed line in a for various By (red: B — 0, scale-bar: 0.1e 2 /h). d-f, The second derivative d 2 I/dV^ d along the center of the 
No + 1, No + 2 and Nq + 3 diamonds, respectively, as a function of a parallel magnetic field. Peaks/dips appear at inflection 
points of the differential conductance and thus correspond to the energy difference between ground and excited states. In f the 
inset shows dl/dVsd vs. —0.3 < V„d < 0.3 mV and By = 0;0.55; 1.1; 1.65 T (arrows) illustrating the splitting and SOI-induced 
reappearance of a zero-bias Kondo resonance, g-i, As d-f but measured as a function of B± . The effective spin-orbit coupling 
appears directly as the avoided crossings indicated by Aso- In d-i the black lines results from the single-particle model with 
parameters Aso = 0.15 meV, Ak k 1 = 0.45 meV, and <? or b = 11.4. The dashed lines in e,h correspond to the excitations to 
the two-electron singlet-like S2 state which cannot be reached by promoting a single electron from the ground state (So) and 
therefore expected to be absent in the measurement. 
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FIG. 4: Tuning Aso in accordance with the curvature- induced spin-orbit splitting of the nanotube Dirac-spectrum. a, Measured 
effective spin-orbit coupling strength as a function of V g extracted from spectroscopy measurements like in Fig. 3, repeated for 
multiple shells. The dashed line is a fit to the theory. Lower inset: Expected dependence of Aso on ejv highlighting the two 
SOFcontributions A so and Ag . b, Graphene dispersion-cones around one if-point of the graphene Brillouin zone. Due to 
SOI the spin-up (blue) and spin-down (red) Dirac cones are split in both the vertical (E) and fex-direction. The cut shows the 
resulting CNT band structure also shown in the upper inset in a. 



